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ON 

°^ 1 Introduction 

The first ideas of a varying Newton's coupling G come from 1937, when P. A. 
M. Dirac introduced his famous Large Numbers Hypothesis [TJ[2]. In the 60s, 
in an attempt to reconciliate Mach Principle with General Relativity (GR), 
Brans and Dicke developed their well known scalar-tensor theory of gravity 
[3]. Following Jordan's ideas, Brans and Dicke generalized GR including a 
varying G, whose dynamics was governed by a scalar field. See [I] for a 
detailed review on varying fundamental constants. 

In addition to the effects of introducing a dynamical coupling constant, 
we know that gravitational interaction is also sensitive to the existence of 
extra dimensions, which could manifest themselves at short distances. In this 
paper, we will be concerned with models that incorporate both a varying G 
and a higher dimensional set up. 

Although the idea of extra dimensions is not new either [SJE], the advent 
of modern (string) theories has brought to the fore the higher dimensional 



1 



scenarios more recently. One of the string inspired models that have attracted 
much attention in the last decade is the Randall-Sundrum model (RS), which 
consists of an effective brane-world embedded in an orbifold of AdS 5 space 
[7]. The property of the RS-like models that is interesting for our purpose is 
the relation between the tension of the brane, and the Newton constant of 
the 4-dimensional effective theory. Models with non-constant brane tension 
thus lead to a time-varying G, as we will discuss below. 

The idea of this paper is to confront particular brane-world models with 
constraints coming form observational cosmological data. The particular 
model we will consider here is that of [8], which is motivated by supergravity 
in singular spaces. We will consider this model as a working example to show 
how observational data could be used to constrain parameters of this type of 
scenarios. 

The observational data we will use to constrain the model are of rather 
different types. For instance, we have data coming from planetary/geological 
scale: Observations due to space missions to Mercury, Mars, Venus and the 
Moon in the 70s, determined that if G varies in time, its variation is less than 
1CT 11 per year. On the other hand, in the late 70s, many works appeared 
relating the relative variation of G with planetary radius [HI [EH]. McElhinny 
et al. studied the evolution of the Earth's radius and extended their study to 
the Moon, Mars and Mercury, and constrained AG/G in specific moments 
close to Solar System formation. 

At cosmological scale, a variation of G leads to modifications in the Fried- 
mann equation. The direct consequences of these variations are changes in 
the cooling rate of the universe and in the computed primordial abundances 
of He and Li. In [TT], Accetta et al. used this relation between varying G 
and light elements abundances to give a bound to the relative variation of 
G. This variation (its absolute value, in fact) turns to be less than 40% since 
Big-Bang nucleosynthesis (BBN). CMB anisotropies are also sensitive to a 
varying G. Chan et al. concluded in [J2] that the relative variation of G 
since recombination is less than 10%. 

In this work, we explore a five-dimensional gravitational model, alia RS, 
with a scalar field in the bulk that modificates the brane tension, which 
induces a variation in G. The variations of G predicted by the model, de- 
pending on two parameters, will be then compared with the observational 
data mentioned above. That is, the aim is to constrain the possible values 
of these parameters, using experimental bounds. 

This work is organized as follows: In Section 2 we discuss the Randall- 
Sundrum-like model coupled to a scalar field, which induces variation of 
effective Newton constant in four dimensions. In Section 3 we survey bounds 
on the variations of G and the observations. In Section 4, we combine obser- 



2 



vational data of Section 3 with the predictions of the model, and use this to 
constrain the parameters. In particular, we give bounds on the 5D-Planck 
mass, supersymmetry breaking scale, and the cosmological constant in the 
bulk. In Section 5, we draw some conclusions. 

2 A brane-world scenario 
2.1 Field equations 

The RS-like scenarios propose that our universe is a 3-brane embedded in 
a curved asymptotically hyperbolic 5-dimensional bulk, or an orbifold of it. 
One can also include matter in the brane as well as in the bulk [S]; here 
we consider a scalar field <fi. The brane is located at the origin of the fifth 
dimension, which we denote x§. This dimension has Zi symmetry in our 
case. 

Fields of the Standard Model live on the brane, while gravitational inter- 
action (and the scalar field) is free to propagate in the bulk. Bulk action is 
given by 

where R is the curvature scalar associated to the 5-dimensional metric g AB ] 
U{4>) is the bulk potential, which depends on the scalar field 0, and 
1 /Mf, being M 5 the Planck mass in 5D. 
The action of the brane depends on its tension Ub(4>) (brane potential). 
In our case, it is a function of the scalar and of the confined matter; namely 

Sbrane = J ^a/"^ 4 ) ( ~ 2k|^ B ^^ 5 = )) + ^matter) , (2) 

with g^s = d^S^g^^^Q. In this paper, Latin indices in capital letters go 
from to 5 (excluding 4), Greek indices go from to 3, and Latin indices in 
regular letters, from 1 to 3. 

The matter content of the 5D space is characterized by the energy- 
momentum tensor, which can be derived from the total action and has the 
bulk and brane contributions; namely 

rp rpbulk I rpbrane fo\ 

AB J- AB ' AB i V 6 ) 

with 

T bu % = I - ^ (5) >(0), (4) 
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T brane \ = (-^ (5) A B UM + r matter ^ 5(* 5 ), (5) 

and 

r matter = diag(-p m , Pm , Pm ,p m ,0). (6) 

Tensor r matter is related to the ordinary matter on the brane. 

The energy density p m and the pressure p m are independent of the position 
in the brane, so one recovers an homogeneous cosmology in four dimensions. 
The equation of state that relates these quantities is taken to be p m = ur m p m . 

Einstein equations reads 

Gab = Rab - ^R 9ab = K l ( T AB k + Tab™) = K lT AB . (7) 
Let us propose the following ansatz for the metric: 

ds 2 = -A 2 (t, x 5 )dt 2 + B 2 (t, x B )dx i dx i + C 2 (t, x 5 )dx 2 5 . (8) 

We are interested in cosmological scale solutions, so we assume an isotropic 
and homogeneous metric in the three spatial coordinates. That is 

ds 2 = a 2 (t, x 5 )b 2 (x 5 )(—dt 2 + dx\) + a 2 (t, x 5 )flijdx l dx^ , (9) 

where flij is the metric of a 3-dimensional space with constant curvature: 

% = ^i(l + jx l x m 5 lm y 2 , (10) 

where the values K = 0, 1, —1 correspond to a (spatially) flat, closed or open 
universe, respectively. Since observational evidences are consistent with a 
spatially flat Universe [14], then we assume K — 0. 

It is important to note that in ^ we chose a conformal gauge for the 
(0—5) part of the metric. In this gauge, the brane is placed in a fixed position, 
£5 = 0, i.e. the fixed point of the Z% symmetry in the fifth dimension. 
Function b only depends on the spatial coordinate x$. 

To derive the brane dynamics, one must verify that, although the equa- 
tions of motion must be restricted to it, these equations have to be satisfied 
in the bulk as well. The brane "proper time" is 

dr = ab\ X5=0 dt, (11) 

and the differential of the normal vector to its surface is given by 

dy = ab\ X5=0 dx 5 . (12) 
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From now on, we write / = /' = J^. 

The Israel-Darmoise junction conditions describe how a brane with a 
given energy-momentum tensor can be embedded in a higher dimensional 
space-time. These equations yield 



y=0 = -4{ P +P), (14) 



V 1 

V y=0 ~ 2 

where equations p is the energy density, and p is the pressure on the brane 
On the other hand, the boundary condition for the scalar is [T5] 



lH 



dU B 

\V=0 Q(f) \y=0 



is) 



The total energy density and pressure on the brane can be written as a 
sum of two contributions: a term related to confined matter, and a second 
one, related to the tension, which in this case depends on the scalar field. 
Thus, we have 

3 3 
P = P™ + ^U B , p = p m - —^U B . (16) 

5 5 

In what follows, all quantities will be evaluated on the brane, i.e. at 
y = 0. Restricting the (0 — 5) component of the Einstein equations to the 



brane, and using boundary conditions (13) and (14) we obtain the energy 
conservation equation 

p = -3H(p + p)-2T% (17) 
where H = -L=o is the Hubble parameter on the brane. 



Using the explicit form for p and p from (16), total energy density con- 
servation law transforms into a conservation law for the energy of the brane; 
that is, 

p m = -3H(p m + p m ). (18) 

Time variation of the scalar field energy density 3£/b/2 cancels the term 
involving T° 5 , since the latter can be written as T° 5 = — \4>'4> = — \U B . 
The solution for the brane energy conservation is 

p m = p a- 3(1+ " m) , (19) 

as in standard cosmology. 
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On the other hand the Friedmann equation on the brane is a consequence 
of the (5 — 5) component of Einstein equations. For a brane containing matter 
coupled to a scalar field 0, Friedmann equation reads 

2 _ K t „2 i K l T T „ 1 f A^.^ Vi2 om K 5 f J_„4„ dJJ-Q A 



21] 



(20) 

with 

The set of equations is completed by the Klein-Gordon equation for the 
scalar field [HHEI; namely 

} + ^ + \(\ ~ ^)P^l = ~% + A*, (22) 

where 

A$- 9V I 9Ub \ 92Ub \ (23) 

Following [16] and [T7], we consider 



A$ = 0. (24) 



Einstein equations have been used to write (22) in this form (see [16j). 



2.2 Physical considerations 



Friedmann equation (20) is not conventional. In contrast to the standard 
one, (20) presents terms that depend on the field 0, a quadratic term in the 
energy density on the brane (present also in absence of the scalar), and an 
additional term that goes like a -4 . 

By the time of primordial nucleosynthesis, corrections coming from brane 
models, including the term proportional to the square of the energy density 
in Friedmann equation, must be negligible. Otherwise, the rate of expansion 
would be modified and the computation of light elements abundances would 
be inconsistent with observations. In this non-conventional scenario, the 
freezing temperature of proton to neutron ratio Tq would be of the order 
of (2 - 3) MeV, while in standard cosmology it is T c ~ (0.7 - 0.8) MeV, 
consistent with He abundance. The difference between both temperatures 
is a direct consequence of the fact that Hubble parameter is linear with T 4 , 
and not with T 2 , generating a slower cooling of the Universe [TS]. However, 
corrections might be important during the inflationary period. 
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Let us be reminded of the fact that in standard cosmology Friedmann 
equation is 

TT 2 S^G A 4 , n „. 

#stand = ^"Pm + y , (25) 

where A 4 is the cosmological constant in four dimensions. Then, the quadratic 



term in p m in (20) can be identified with the first term in (25); that is 



tm 4 = **g 

12 5 3 y ' 

It is clear that in our model, Newton constant in 4D varies as it depends 
on 0; i.e. it is possible to find time variation of G. 

2.2.1 Bulk and brane potentials 

We consider a functional form for the potential U((f>) coming from the su- 
pergravity models in singular spaces studied in [15]. Following these results, 
one finds 



u = (s^J - w ■ (27 » 

where W(4>) is the so called superpotential. 

We study the case in which the superpotential in an exponential function 
of the field 

W{<j>) = 4:ke a *, (28) 

where [A; -1 ] = L and a is a real number. 

The brane potential is defined through the superpotential by 

U B = TW } (29) 

where T is a real number related to the scale of supersymmetry breaking 

Having the functional relation between U B and the scalar field, given by 



(28) and (29), one can find G{4>) using (26); thus, 

G(0) = £-%Te°*. (30) 

The expectation value of today is assumed to be zero by convention as 
a boundary condition. Then, Newton "constant" would be given by 

Gtoday(0) = £-*%r. (31) 
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2.2.2 Working hypothesis 

The form of the Friedmann equation with all the new contributions is quite 
abstruse. Then, in order to solve the model, some approximations and as- 
sumptions have to be taken into account. We discuss these below. 

First, the term proportional to a -4 , can be considered as a correction to 
the radiation density. This term is usually referred as dark radiation. Here, 



we assume A = in (20). We also consider a low energy regime, i.e. we 
neglect the term proportional to in (20). It is possible to do this under 
the condition p m -C p C Ht, with 

Pc ri t = ^ = ^T~4.6xl0 33 ^, (32) 
k\ ki cm 6 

.2 



where (28), (29), and bounds on k, k\ and T consistently found a posteriori 
in Sections 3 and 4, have been used. In the studied period (between BBN 
and today) the density remains below this critical density. We also assume 
that the time evolution of the scalar field in the brane proper time r is 
much slower than the one of the scale factor a. It is possible to extract 4> 2 



from the first integral in (20) in this adiabatic regime. 

A non-dissipative approximation of the potential will be also considered. 
The brane potential U-q is basically Newton constant on the brane, up to 
multiplicative constants. Following the adiabatic approximation, it is rea- 
sonable to suppose that the contribution correspondent to this term might 
be negligible. The term dUs/dr, as well as the other terms of order <p and a, 
contribute to higher order estimation of G(<f>) than the one we study here. 

Finally, and consequently with the assumptions above, the square of the 
time derivative of 0, i.e. the kinetic energy of the field, is lower than other 
terms in Friedmann equation. It is possible to make a simply calculation to 
constrain the current value of the time derivative of the scalar field: Following 
the approximations, Friedmann equation today, divided by Hq is 

i = n M + n R + n A + (33) 

where VLq^ = 2 (r o )/16. 

Recent data [O] implies that the sum of the first three contributions is 
close to 0.996, which fixes a limit to the absolute current value of the time 
derivative of 0; namely 

|A( To )| < .22. (34) 

This heuristic argument supports our approximation hypothesis. 

In addition, the exponential dependences appearing in Friedmann and 
Klein-Gordon equations will be approximated to 1 since <fi is small. 
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With all the approximations described above, Friedmann equation takes 
the form 

H 2 = T k -4 Pm + (35) 

with 

^ :=k 2 (T 2 -l)(l-a 2 ). (36) 



On the other hand, eq.(19) says that, for radiation, 

PR = p R#o a ~ 4 > ( 37 ) 

and for non-relativistic matter, 

Pm = PomHqCl~ 3 , (38) 

with Pom,r dimensionless constants. 

Thus, Friedmann equation can be rewritten as follows: 

H2 = T r*(^ + P J) HA ° + ir- (39) 

After making the identifications 

k 

T-kIpom,rHI = fi M ,R, (40) 

^r n - < 41 > 

we obtain the familiar form for the equation 

H 2 „ / \ 

wr^ + ^ +a "- (42) 

Then, the equation for <fi is 

4> + AH^ + 2(^-u m )p m 4akT=-16a^. (43) 

The system of equations above can now be solved in two different epochs: 
one dominated by radiation and matter, and the other governed by matter 
and cosmological constant. 
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2.3 Solution to the Field equations 
2,3.1 Radiation and matter dominated epoch 

In this epoch the term proportional to A e g can be neglected. Then, the 
Fridmann equation takes the form[j] 



whose solution is 



a i(v) = ^(V-Vo) 2 ~ (45) 

where the integration has been performed in the conformal time r], being 

dr= a -jg>d V . 

The equation of state for radiation is p ra d = Prad/3, while p mat = is the 
equation for non-relativistic dust. Then, KG equation reads: 

4> x + 4#i0i + 2aHl^- = -l6aH$tt A , (46) 



a 



where we used eqs. (40) y (41). The solution for the differential equation is 
(see Appendix) 

M"i)=B-—^>-^"i, ( 4 7) 

where it has been taken into account that scale factor remains small during 
this regime. 

2,3.2 Matter and cosmological constant dominated epoch 

During this regime, the term proportional to the inverse of the fourth power 
of the scale factor a is negligible when compared to the other two terms. 
Then, the Friedmann equation reac 1 ^' 



§| = ^ + n A . (48) 

Integrating, one obtains the scale factor during this epoch, namely 

a 2 (r) = a sinh§ ( ^H (r - C )) . (49) 



1 Subindex 1 refers to the radiation- matter regime. 
2 Subindex 2 refers to the matter- A regime. 
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In this regime, the Klein-Gordon equation is 

2 ^M 



+ 4H 2 (f) + 2aH'f-^- = -lQaHfo A . (50) 
a 2 



In order to solve this equation it is convenient to consider two cases: 
a 2 <C 1, that occurs for times close to union time (r ~ r u ); and a 2 ~ 1, it is, 
times near today (r ~ r ): 



a 2 < 1: 



^2(02) = D - 



2C 



5\Z^M^2 



- -Qtln(o 2 ; 
5 



(511 



02 ~ 1: 



i(a 2 ) = G 



{8n A + n M )al, (52) 



where G and F are integration constants. Details of this calculation can be 
found in the Appendix. 

2,3.3 Boundary conditions for the scale factor 



To find the integration constants in (45) and (49), one must define the bound 



ary conditions. We take the convention 

02(^0) = 1, (53) 

where To is the current value of r. We also have the matching conditions 

H 2 (r ) = H , (54) 

and 

#i(r tt ) = H 2 (r a ), (55) 

where r u is the junction time between both regimes. In addition, we have 

ai(r n ) = a 2 (r u ) = a u . (56) 



Constants a and C can be obtained from (53) and (54). Consequently, 
we have 



a 2 (T) 



1 1 V • , 2 / 3^A rr / 2 . , _j 
— 1 smhs if t — t H ;= smh 



l-^A 

(57) 
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The constant of integration must be chosen to satisfy the condition 
ai (r = 0) = 0. Then 



T-(di) 



3 

2 4^ 2 

(Q M Qi ~ 2^ R ) a/ ^m«i + + o o- o2 • 



(58) 




t[10 years] 



0.01 



0.05 



0.09 




Figure 1: Graphics for r vs ai and r vs a 2 



2,3.4 Boundary conditions for the scalar field 

Since the scalar field must be smooth, (pi (and its derivative) must be equal 
to 2 (resp. to its derivative) in the junction time r u , when 
a = a a = flji/fl\. That is, 

0i(a u ) = 2 (a u ), (59) 

^K) = ^K). (60) 

The field 2 must be equal to 3 (as well as their derivatives) in an 
intermediate time T\ between r u and To. We take T\ as the time equidistant 
to r u and r . Thus, 

2 (ai) = 03(ai), (61) 

^ = (62) 

The last boundary conditions we need are the values of the field and its 
derivative today, i.e. at r (o(t ) = ao — 1) : 

0s(ao)=O, (63) 
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da>2 



po- 



rn 



Using (59)-(64) and the values for H , fi M , and f2 A listed in Subsection 
3 of the Appendix, one obtains the constants of integration A, B, C, D, F and 
G. Table Q] shows the values of the constants as linear combinations of a and 
Po- 



Constant 


a 


Po 


A 


0.09469 


0.0309 


B 


682.718 


57.114 


C 


1.745 


0.719 


D 


2.212 


0.752 


F 


2.408 


0.998 


G 


2.0104 


0.333 



Table 1: Values of the integration constants of the fields which are linear 
combinations of the parameters a and po. 



The final expression for the fields can be found in the Appendix. In 
Figure [|J the behavior of the field for different allowed values^ of a and p 
are observed. Figure ^ A) corresponds to a = 1/V20078 and p = —0.02164, 
which is the value for p when a has the pointed value. Figure shows 
the graphic for the field, for a null value of a, and po = 0.22, which is the 
upper limit for po coming from Friedmann equation. Finally, Figure [2|C) 
shows the unique case in which the field does not diverge at the origin, due 



to the fact that the value of a is such that the linear term in a in 



eq.(102) in the Appendix) is zero for all po- 



>i (see 





<* (C) 



-10 



-20 



-40 




Figure 2: Different behavior of <\> vs a. (A) a = 1 / v / 20078 and p = -0.02164; 
(B) a = y po = 0.22; (C) a = -0.326062p and p = 0.22. In this case the 
field has no divergence at the origin, and its value there is -36.4087. 



3 Allowed values for a y po are detailed in Section 4. 
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3 Experimental and observational data on G((f>) 



3.1 Corrections to Newtonian potential 

Theoretical speculations predict new effects at distances of order less than 
1mm. In particular, models with spatial non-compact extra dimensions are 
of interest because these "internal" dimensions could alter the form of the 
Newtonian potential. 

If the extra dimension is non-compact, as in the case of RS model, there 
is a continuous of Kaluza-Klein (KK) modes for the gravitational field. The 
continuous spectra of KK modes leads to a correction to the force between two 
static masses in the brane. The potential for two point-like masses confined 
to the brane reads [7] 

V5a(r) = G=!^fl + JL\ <«) 



where k' 1 should be of order of the distance of available gravitational tests 
(~lmm ) or smaller. 

Thus, in order to bound the deviation from the Newtonian potential, one 
must constrain parameter k. Adelberger et al. [2U] performed experiments 
with torsion balances to model the correction to Newtonian potential using 
a power law of the form 



(66) 



where the values of j and \(3j\ are shown in Table [2 



k 


\m<) 


2 


4.5 x 10" 


-4 


3 


1.3 x 10" 


-4 


4 


4.9 x 10" 


-5 


5 


1.5 x 10" 


-5 



Table 2: Bounds on \/3j\ for j = 2,3,4,5 obtained by Adelberger et al. [20j. 



For j = 3 the extra term is 
while the deviation predicted by RS model is 

,77117712 f 1 



at/RS 



G- 



r 2 k 2 



(67) 



(68) 
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Thus, the value of k is constrained comparing AVil with AV^ S . Then, 
the characteristic scale at which the effects due to the presence of an extra 
dimension become important is 

- < 0.01mm. (69) 

3.2 Observational bounds on G variation 

Bounds on the variation of the Newton constant are obtained from local 
and cosmological observations. Local observations are related to the solar 
system, as well as nearby stars. Geological and paleontological data, as well 
as planetary orbits, stellar densities, and luminosities, are of great importance 
when studying G, because they are affected by its variation. 



3.2.1 Bounds on AG/G 

Planetary radius variation: In 1961, Egyed proposed that paleomagnetic 
data could be used for the calculation of Earth paleoradius (past to current 
planetary radius ratio) in different geological eras. Starting from the hypoth- 
esis that the continental material area remained constant during planetary 
expansion, Egyed found that the ratio between current and past angular sep- 
aration (paleolatitud) of two given sites is proportional to the paleoradius [9j. 
A few years later, in 1978, McElhinny el al. related Earth radius variation 
to time evolution of gravitational constant, and extended the analysis to the 
Moon, Mars and Mercury |T0]. According to their work, 

AR AG 

~R~ = ~"'~G~' (70) 

where AR is the variation of the radius R and 7 is a constant that depends 
on the planet structure. 

On the other hand, there is another way to write this relation using the 
paleoradius i? a : 

Cx 7 

Table [3] summarizes the results for the Earth, the Moon, Mara^and Mer- 
cury. These results assumes that the surface of each studied planet acquired 
its current shape by the time indicated in the fourth column. 



4 There are two different analysis for Mars: (A) assumes a 19km expansion during the 
last 3600 million years; (B) supposes a 1km variation on martian radius in the past 1000 
million years. See [10] for details. 
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Planet 


R a 


7 


Time [10 9 years] 


|AG/G|(<) 


Earth 


1.020 ±0.028 


0.085 ±0.02 


0.4 


0.62 


Moon 


1.0000 ±0.0006 


0.0004 ±0.001 


3.9 


1.5 


Mars (A) 


0.9944 


0.03 ±0.01 


3.6 


0.12 


Mars (B) 


1.0000 ±0.0003 


0.03 ±0.01 


1.0 


0.01 


Mercury 


1.0000 ± 0.0004 


0.02 ±0.005 


3.5 


0.02 



Table 3: Paleoradius (R a ) of the Earth, the Moon, Mars (A and B) and 
Mercury, and bounds on relative variation of G. 



Big-Bang Nucleosynthesis: Bounds of a different sort come from cosmol- 
ogy. In 1990, Accetta et al. studied bounds on gravitational constant value 
during primordial nucleosynthesis, considering neutron mean life measure- 
ments [H]. They determined D, 3 He and 7 Li abundances while varying G, 
and how this variation affects barion to photon ratio. On the other hand, 
Copi et al. recalculated relative variation of G since BBN, but using only 
primordial D abundance in quasars [2T] . In both works the constraint on 
relative variation of G is 

AG 

|-^|bbn<0.4, (72) 

which means that the relative variation of gravitational constant since BBN 
is less than 40%, at the 95% confidence leveQ 

Cosmic background anisotropics: Power spectra of cosmic microwave 
background anisotropies (CMBA) can be useful while constraining G vari- 
ations in cosmological scales. In [T2] . gravitational constant stabilization 
(convergence to its current value) and its relation to CMBA are studied in 
detail. Two possible parametrizations of G are considered: one, correspond- 
ing to an instantaneous stabilization, and the other to a stabilization linear 
with the scale factor a. 

If the stabilization is linear, the relative variation of G since recombination 
(z ~ 1000) is 

|-^-|cmb<0.1. (73) 
at the 95% confidence leveQ Thus, the relative variation of G (its absolute 



5 It is important to note that the constraint is for the absolute value of ^^Ibbn- 
In [21] the constraints are: 0.85 < Spt < 1.21, at the 68.3% confidence level, and 

0.71 < ^q^- < 1-43, at the 95% confidence level. G is the present value of the Newton 
constant. 

In this work we made use of the last constraint. 

6 Again, we aware the reader that the constraint is for the absolute value of ^P|cmb- 
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value) since recombination is less than 10%. 



3.2.2 Bounds on G/G 

Lunar Laser Ranging (LLR) has been measuring the position of the Moon 
with respect to the Earth during more than thirty years, with a precision 
of 1cm. The misions Appolo 11, 14 y 15, and russian-french Lunokhod 1 y 
4 carried retro-reflectors to the Moon, which reflect laser pulses sent from 
the Earth. LLR data are used to constrain Weak Equivalence Principle, 
post-newtonian parameters and G/G. 

According to 2004 data in [22], the maximum variation allowed to gravi- 
tational constant today is 

^(r ) = (4±9)xl0- 13 yr- 1 . (74) 



3.2.3 Bound on G/G 

Now, let us discuss the constraints coming from observational bounds on 
G/G. For that purpose, we consider a model slightly different to brane 
cosmology. 

Other model which considers scalar fields is the scalar-tensorial theory of 
Brans and Dicke of 1961. This theory contains a scalar field governing G 
dynamics. In usual notation we have |3], 

Gfo>) = -. (75) 

G time dependence is described by 

G(t) ~ r" n , (76) 

where n = 2/(4 + 3u>), and u is a model parameteiQ which measures the 
deviation from General Relativity (GR). GR results are reobtained when u> 
goes to infinity. Considering ([761, it can be shown that 



|(r) = n(n + l)r- 2 . (77) 



In [12] the constraints are: 0.95 < G ™ B < 1.05, for a G variation modeled by a step 
function, and 0.89 < G °*m <- 1.13, for a variation modeled by a linear function of the 
scale factor. Both constraints are at the 95% confidence level. 
In this work we made use of the last constraint. 

7 In Brans-Dicke work ui = const., but there are more complex models where u> = oj((p). 
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Benvenuto, Althaus and Torres in [23j gave a bound to the absolute value 
of u, coming from white dwarfs evolution and its relation with a varying G. 
According to their results, the calculated luminosities differ from the observed 
ones for |u;| < 5000. Then, the allowed values range is |u;| > 5000. 



Equation (77) evaluated today (r =13730 million years) together with the 
constrain on u, give a bound on the variation of the second time derivative 
of G: 

G 

- 3.55648 x 10~ 4 V 2 < — (t ) < 7.11082 x 10~ 4 V 2 . (78) 

Gr 

At this point, one could ask about the relation between Brans-Dicke 
model and the one studied in this work. The explanation is the following: 
Brans-Dicke theory in the limit us — > oo gives back GR; while GR and brane 
cosmology should be equivalent in the studied limit. Then, at first order, BD 
and brane cosmology are equivalent and constraints on oj can be translated 
into constrains on brane model parameters. However, it is not clear that the 
differences between BD and brane cosmology lead to a great discrepancy on 



the limit (78). Then, we study the inclusion and exclusion of this bound in 



following analysis. 

4 Constraining model parameters 
4.1 Parameters a and 0(to) 

G variations, i.e. AG/G, G/G and G/G, can be written in terms of as 
follows: 

a) = —ac(f)(a), (79) 
a) = a0(a), (80) 



AG 
~G~ 

G 
G 



(a) = a 2 4> 2 (a) + a<f>(a). (81) 



G 
G { 

Then, observational bounds discussed above restrict the possible values of 
a and pq. The strongest restrictions on a and po come from BBN and recom- 
bination. The largest value for the modulus of a is fixed by the combination 
of BBN and CMB restrictions; that is 



M < ; = 0.007, (82) 

v / 20078 



while \p \ < 0.22. 
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The allowed values are shown in Figure [3] Allowed values are inside the 
surface which perimeter is given by the green hyperboles. 

Po Po 

0.4 

0.2 
o 

-0.2 



-0.006-0.004-0.002 0.002 0.004 0.006 -4x10 _9 -2x10 9 2x10 9 4X10 -9 

Figure 3: Allowed values for a and po. The figure on the right is a zoom- 
in of the left one. Allowed values are located inside the dark gray contour, 
and satisfy |a| < 0.007 y \p \ < 0.22. Blue lines in the left figure establish 
the bound for the largest value for \p \, a condition coming from Friedmann 
equation today. 



4x10" 8 
2x10 8 

o 

-2x10 8 
-4x10 8 




In [21], Brax and Davis find two theoretical values for a, emerging from 
supergravity in singular spaces. This values are a = l/\/3 and a = — 1/VT2. 
From the analysis we have just done, based on observational constraints, the 



absolute value of a is bounded by (82) and then, it excludes these theoretical 
values. 

4.1.1 Statistical analysis 

Now, we would like to investigate how robust our constraints are. That is, 
we will study how the constraints change if one excludes ones or others data. 

Excluding G relative variation since Big-Bang nucleosynthesis data, and 
using the other data, i.e. AG/G| C mb, AG/Glp^^a, G/G\ todSi y, G/G\ today 
and \p \ < 0.22, it can be deduced that the strongest constraint on a absolute 
value is given by the combination of AG/G\qmb and G/G\t da.y 

In this case, we have 

H < -= = 0.024. (83) 

v / T787 

If we exclude CMB data, in addition to BBN data, and make a similar 
analysis, we have that the constraint on \ol\ is fixed by the combination of 
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G/Gjtoday and G/G|today: 



\a\ < 



V237 



0.065. 



(84) 



Omitting observational data of the second derivative of G, we obtain the 
same results that those at the beginning of this subsection. This is because 
the severest constraint on |a| is given by the combination of BBN and CMB 
data. 

If we exclude BBN data, in addition to G/G|today ? a is constrained by the 
combination of AG/G\cmb and G/Gjtoday cotes, being 

' (85) 



a < 



V243 



0.064. 



Finally, if we only take into account the bounds coming from today's 
value of G derivative, and those related to planetary paleoradius (i.e. BBN, 
CMB and G/G\ today limits, excluded), \a\ has no bound. 

Table [4] condenses the results obtained in data analysis. 



\a\ 


Obtained with: 


l/y/3 = 0.577 
1/VT2 = 0.289 
1/V237 = 0.065 
1/V243 = 0.064 
1/V1787 = 0.024 
1/^20078 = 0.007 
no limit 


SUGRA 
SUGRA 

G/G today; G/G\ today 

CMB; GyG|t oday 
CMB; G/G|today 
BBN; CMB 
G/G today; Mercury paleoradius 



Table 4: Constraints on a obtained combining G variations observational 
data. The first two values are equalities found with supergravity (SUGRA) 
in singular spaces. The other values are cotes and should be read "|a| < 
The severest constraint is given by Big-Bang nucleosynthesis (BBN) and 
CMB. All the constrains, except the last one, exclude the first two values. 



In the studied cases, the upper bound on \p \ is fixed by Friedmann equa- 
tion evaluated today, i.e. |p | < 0.22. In addition, it should be said that |a| 
is always below 0.065, except in the last case, where its value has no bound. 

Analog analysis can be done, but assuming a different bound on BD 
parameter: \uo\ < 500 (instead of \oj\ < 5000). Table [5] shows the results. 

The difference with respect to the results obtained with |u;| < 5000 lies 
in the constraint of |a| found combining CMB data and G/G|t day bound. In 
this case, the constraint on |a| is less restrictive since it is three times bigger. 
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a 




Obtained with 


1/V231 = 0.066 
1/^243 = 0.064 
1/V178 = 0.075 
1/^20078 = 0.007 
no limit 


G/G today; G/G today 

CMB; G/G|today 
CMB; G/G|today 
BBN; CMB 
G/G today; Mercury paleoradius 



Table 5: Constrains on a when combining observational data of G variation, 
in the case M < 500. 



4.2 5D parameters, T and A e // 

The 5-dimensional parameters of the model can also be constrained, as well 
as the 4-D cosmological constant A eff . Using H and from Subsection 3 of 



the Appendix, the constraint on k from (69), and |a| < l/v20078, we have 
the results of Table [6] 



Parameter 


Bound 


«? 


< 2.8 x 10~"s 3 


M 5 = 1/4 


> 4.7 x 10 s GeV 


As 


< -4.2 x 10 2 V 2 


\T\ 


< 1 + 3.05 x 10" 63 


A eff 


= 4.8 x 10" 18 s- 2 



Table 6: Constrains on 5D parameters, T and A e ff. 



The same results are obtained when the other constraints on a from 
Tables [I] and [5] are used. Constraints on k 2 and M 5 are in accord with the 
ones predicted in [251 1261 127] . 



5 Conclusions 

In this work, we studied a brane-world cosmological model in which varia- 
tion of the Newton coupling G emerges naturally. This model is inspired in 
supergravity in singular spaces [21] ■ By resorting to available observational 
data, we manage to constrain the parameters of the theory. Light elements 



abundances, coming from Big Bang nucleosynthesis (He, Li y D)(eq. (72)) 



CMB (eq. (73)) and, near in time, planetary radii variations, allowed to con- 
strain the relative variation of G, i.e. AG/G (see Table [3j. Measurements 
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of Lunar position with respect to the Earth with LLR offered us a bound 



on today's value of the G time derivative to G ratio (eq. (74)). Combining 
the severest constraints on G variations, i.e. those coming from Big Bang 
nucleosynthesis and CMB, a bound for the absolute value of the parameter a 
was obtained. In fact, this parameter must be less than 0.007, and the value 
of |po| is bounded while evaluating Friedmann equation today. 

Statistical analysis was performed on the results to analyze how robust 
the bounds are against the exclusion of particular set of data. That is, we 
studied how the upper bound of \a\ gets affected if one excludes different 
sets of data. Results are shown in Tables [4] and [5j It is worth mentioning 
that the upper bound for |a| is always less than 1/V178 ~ 0.07 (except in 
the case where Big Bang nucleosynthesis and CMB data are excluded, which 
turn out to be the most important ones). 

Our analysis presents a method to investigate the phenomenological vi- 
ability of models that, among other features, predict time variation of the 
fundamental couplings. It could be interesting to extend our analysis to other 
brane-world type scenarios. 
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A Field equations for the scalar field 

A.l Matter and radiation regime 

During matter and radiation epoch. Klein-Gordon equation reads 

0i + 4#i0i + 2aHf^ = -lQaHfo A , 



where we used (40) and (41) 



Differentiating with respect to a\, we have 

#i , 4 -2aH ftt u , 
h —-01 = , == [—2 + 8Qa q i , (87) 



where we used (44) and ipi — <pi- 

Being reminded of the fact that during this regime ai is small, the func- 
tion on the right hand side can be approximated by its first term in power 
expansion for a\ ~ 0, thus 
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1 fa = -2aH - 



n 



■M 



Integrating, we find 



dai ai V^Rti 



■2 ' 



Wi) = a fcT + ^^' ( 89 ) 

where A is an integration constant. 

To find 0i (ax), it is necessary to integrate once again, taking into account 



that 

d(j>i 1 d0i 1 



Then, we have 

0i( a i) = B + / daiip(ai) Ql = ~ £> + / daiip(ai) — Q * , (91) 

J ifoV^MQl + J H y/Q R 

Where B is a constant, and where we only considered the first term in the 
power expansion of "TnH-^^n^ - 

The solution for the differential equation is 

M"i) = B - —*=— - 2 ^ ai . (92) 

A. 2 Matter and cosmological constant regime 

During this regime the Klein-Gordon equation is 

<f> + 4H 2 <j) + 2a# 2 % = -lQaHl fi A (93) 
a 2 

which, written function of a 2 derivatives, is 

di) 4 , —2aH _ ^ 9 , , 
-^ + -^=- 9 (fiM + 8fi A a 2 ). (94) 

dfl2 ° 2 a 2 VO M + ^Aa| 



First, consider a 2 C 1: 



(fi M + 8fi A a 2 )-^. (95) 

_ 9 



«2V^M + ^Aa| a 

Replacing this equation in (94), we have 
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ip2{a 2 ) = 3 + — |-, (96) 



5a| a 2 



where C is an integration constant. 

The solution for the field comes form the integral 

2 (a 2 ) = D+ [ da 2 ip 2 (a 2 )—^— ^ D + / da 2 ij) 2 (02)77" A /7r~> ( 97 ) 

y a 2 -«2 y -no v "m 

where D is a constant, and (48) and (96) were used. Also, l/a 2 H 2 was 
approximated by its first order in the power expansion, for a 2 <C 1. 
After integrating, we have 

2C 4 

2 (a 2 ) = D j - -a ln(a 2 ). (98) 

Now, consider a 2 ~ 1: 

, ; 1 (Am + 8^) * 8 ^±^ . (99) 



As above, and making the approximation 

1 1 



a 2 H 2 H y/Q\ + Q M 
we obtain the approximate value for the field 



(100) 



where G and F are constants. 



A. 3 Final form for the field 

The values we take for the cosmological parameters appearing in the ex- 
pression of the fields are H = 22.69 x lO^V 1 ; fi M = 0.28; fi A = 0.716; 
Q R = 4.6 x 10~ 5 (see [H] for further details). 

The values for the integration constants can be found with the parameters 
above and with the boundary conditions discussed in Section 2. The value 
for the scale factor in the intermediate time a\ is 0.542. Then, 
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, s 6.9539a + 2.2674p , v 

' h (ai) = 682.718a + 57.1142p ~ 2 — - 4029.68aai, (102) 



2 (a 2 ) = -(543.028a + 224.698p ) + (15241.3a + 6283.55p )(a2 - a u )- 

- (297480a + 122591p )(a2 - a a ) 2 , (103) 

03(^2) = Po(a 2 - a ) - (6.0313a + 2p )(a 2 - a ) 2 , (104) 

with a u = 0.0897 and a = 1. 

Notice that 2 series is around a u = 0.0897, while 3 series is around 
a = I. 
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